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Abstract 

We investigate whether postulating the existence of Lorentz-violating, CPT-conserving interac- 
tions allows three-neutrino solutions to the LSND anomaly that are also consistent with all other 
neutrino data. We show that Lorentz-violating interactions that couple only to one of the active 
neutrinos have the right properties to explain all the data. The details of the data make this 
solution unattractive. We find, for example, that a highly non-trivial energy dependence of the 
Lorentz-violating interactions is required. 



I. INTRODUCTION 



In 1996, the LSND collaboration first presented statistically significant evidence for v^. 



appearance from antimuon decay at rest 



Further analysis of larger LSND data samples 



(which also included data from pion decay in flight) confirmed the original excess Q]. In spite 
of the fact that this so-called LSND anomaly is close to completing ten years of existence, 
there is no compelling physics explanation for it. It is very natural to try to understand the 
LSND anomaly by postulating that there is a probability for what was originally a muon 
antineutrino to be detected as an electron antineutrino. If this is the explanation for the 



LSND anomaly, the associated transition probability need not be very large: P„e ~ 0.2% 



Neutrino oscillations are the leading candidate explanation for the LSND anomaly, for 
several reasons. First, neutrino oscillations are known to occur in nature - it is all but 
certain that mass- induced neutrino oscillations are the solution to the solar and atmospheric 
neutrino puzzles Q]. Second, the L dependence of the — > v^, oscillation probability 
is crucial for reconciling the LSND anomaly with negative results from the KARMEN-2 
experiment j^. 

Mass-induced oscillations among the three known "active" neutrino species, however, are 
known not to provide a viable solution to the LSND anomaly. The reason is very easy to 
understand. In terms of two-flavor z/^ v^. oscillations, the LSND data (combined with 
other "short baseline" data) require Am^ ~ 1 eV^ and sin^ 2^e/^ ~ 10^^. On the other hand, 
the well-established atmospheric and solar neutrino oscillations require Am^ ~ 10^'^ eV^ and 
Am^ ~ 10~^ eV^, respectively. With three neutrino masses, it is impossible to obtain three 
mass-squared differences that differ by orders of magnitude. For this reason, the existence 
of so-called light sterile neutrinos (that couple to the active neutrinos via mixing) is often 
postulated when it comes to addressing the LSND anomaly. It is remarkable, however, 
that the addition of one sterile neutrino does not provide a good flt to all neutrino data 
fl a, wKHe the add,«o. of two (o. .no.) ste.le .e„tnnos .ee,n. to p.ovde a bette. fit 0. 
The goodness-of-flt of 3 -|- 2 models when one includes all of the world's neutrino data is, 
however, currently under detailed investigation,^ and, furthermore, concordance cosmology 
fits to all available "cosmological data" virtually rule out sterile neutrinos masses larger than 



^ We thank Michel Sorel and Osamu Yasuda for bringing out this point. 
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^ 1 eV if these mix significantly with the active neutrinos. This is precisely the case of 
sterile neutrinos that help resolve the LSND anomaly 

The fact that mass-induced oscillations do not provide a compelling solution to the LSND 
anomaly is best illustrate d by the fact that there are several alternative solutions to the 
LSND anomaly S H Q Q- While many of them [8| are now ruled out by more 
current data or more refined data analysis (see for example, 0, Q, 15 1), most of the ones 
that seem to work require new, very light degrees of freedom (say, sterile neutrinos or 
scalars) and extraordinary new physics (say, CPT violation) . Potential exceptions include 
postulating that CPT-invariance is violated and that neutrino propagation is not described 



by the Schrodinger equation |lO 
may direction-dependent effects 



Lor irnposing Lorentz-violating interactions that include 
111 . Il^ lisj l . It is yet to be properly demonstrated whether 
either of these two latter possibilities can indeed accommodate all the neutrino data. 

It is interesting to explore whether one could resolve - at least in principle - the LSND 
anomaly without introducing any new, light degrees of freedom. Since new parameters need 
to be added to the ones present in the z/SM,^ we will postulate that one of the neutrinos cou- 
ples to a source of Lorentz Invariance Violation (LIV). By varying the flavor composition of 
this neutrino and the energy dependence of the LIV effect, we are able to probe for solutions 
that qualitatively accommodate all the neutrino data. The reason for this particular choice 
of new physics is straightforward, and will become clear in the next sections. It provides a 
minimal number of free parameters to fit the LSND data, and also contains a mechanism 
for allowing the rest of the neutrino data to be properly interpreted in terms of "ordinary" 
mass-induced neutrino oscillations. 

Our intentions are two-fold. First, as already stated above, we wish to determine whether 
there is indeed such a solution, and how constrained it is. We found that when only one 
superposition of neutrinos is coupled to the Lorentz violating sector the general characteristic 
of the data can be reproduced. That is, the Solar neutrino and Atmospheric neutrino are 
describe to a very good approximation by the z/SM. LSND data, however, is taken care of 
by the new LIV interaction. Second, we will use our LIV solution to estimate how much 
"fine-tuning" is required for a three-flavor solutions to the LSND anomaly. Here, our main 



^ The i/SM refers to the standard model of particle physics, plus the addition of mass parameters for the 
Ught neutrinos Q|. 
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result is rather negative. We find that in order to accommodate all neutrino data, the LIV 
model has to be highly finely-tuned. Thus, while we do report here on our attempts to find 
a LIV solution to all neutrino data, our conclusion is that such a solution is unlikely to be 
realized in Nature. 

The paper is organized as follows. In SeclHl we present the mechanism we wish to explore, 
and how it modifies neutrino oscillations. In Sec. IIIH we qualitatively fit all neutrino data 
with the model outhned in Sec. HIl In Sec. IIVI we speculate on possible origins for LIV and 
summarize our results. 

Our study is qualitatively different compared to other studies of Lorentz invariance vio- 

nnn n n 

lation and neutrino oscillation experiments [Ul, llj, ll3l • R-sfs. jll[ and [13[ concentrate on 
observable consequences of LIV for "short" baseline neutrino experiments, including LSND. 
(Ref. very briefly comments on other neutrino oscillation experiments.) In Ref. [l^, an 
attempt is made to flt all neutrino data. Various different scenarios of Lorentz invariance vi- 
olation are discussed, but the flt is done for subsets of the neutrino data at a time. Moreover, 
only oscillation with L/ L, and LE dependent were considered. Based on our results, we 
believe that none of the examples explored in can flt all neutrino data once a global 
flt is performed. This conclusion is based on the energy dependency, regardless of whether 
neutrinos are assumed to be massive, whether the LIV effects are directionally dependent, 
or whether CPT is conserved. In contrast, our emphasis is on trying to qualitatively ac- 
commodate all neutrino data, including those from LSND. Unlike other studies, our main 
idea is to couple only one neutrino to the LIV source. As a specific model we concentrate 
on CPT-conserving, direction-independent effects. To fit all data we find ourselves forced 
to explore more exotic energy dependent LIV effects. 



II. FORMALISM 



We postulate that one of the neutrinos couples to a source of Lorentz-invariance violation 
(LIV). This assumption is motivated by minimality but turns out to be very important, as 
will be discussed later. 

In the case of one neutrino with mass m, we postulate that LIV manifests itself via a 
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modified dispersion relation of tlie form 



m 



E r^\p\ + — + 



film 



(1) 



2\p\ 2\p\ 

Eq. (0) is valid as long as \p\ ^ m, / is some function of so that CPT is conserved.^ 
We assume that / can be expanded as 



2n 



(2) 



where Eq is some convenient constant with dimensions of energy, while the coefficients 6„ 
have dimensions of mass-squared. 

Assuming that all approximations above are valid, neutrino oscillations are described by 
the effective Schrodinger equation 



. d 

dL ° 



(3) 



where a, P = e, fi,T and 
H = U Ai2 



n=l 



cosC cos6'i^^ 



COS C sin 6l 
\ sin C ) 



cos C cos 6*^ cos C sin 6^^ sin^^- 



(4) 

Here, A,j = Amjj/2E, i,j = 1,2,3 where Am]j are the neutrino mass-squared differences 
and U is the leptonic mixing matrix, while a„ = hn/ Eq (the a„ coefficients have dimensions 
of mass). The mixing angles C, and 9^ are defined by Eq. (^. They characterize the flavor 
content of vl, the neutrino that couples to the LIV sector. We have assumed, for simplicity, 
that there are no CP-odd parameters in the LIV sector. 

All information we need to analyze the neutrino data is included in Eq. (j^. Before 
proceeding, it is convenient to describe how we hope to accommodate all the neutrino data: 

• For large energies and short distances, |Ajj|L <^ 1, all oscillation behavior is governed 
by the LIV terms. These will be flxed by the LSND anomaly plus other "short- 
baseline" constraints. 



^ More generally, we could have chosen / to also contain odd powers of These, however, violate CPT, 
and lead to a different modified dispersion relation for neutrinos and antineutrinos. We find that this extra 
complication does not add significantly to our analysis, and we choose to work only with even powers of 
\p\. This choice is natural in the sense that it is protected by a symmetry (CPT-invariance). 
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• For "solar" oscillations, we will take advantage of the energy dependence of the LIV 
term in order to make sure that these are negligible at typical solar neutrino energies. 
In particular, we will find that all a„ will vanish for sufficiently small values of n. 

• For "atmospheric" oscillations, LIV effects will turn out to be very small. At large 
enough energies, the neutrino that couples to the LIV sector is a Hamiltonian eigen- 
state, while, because the LIV part of the Hamiltonian has rank one (two zero eigen- 
values), the other two eigenvalues arc still of order A12, A13. We will choose ( and 6^ 
so that the "LIV neutrino" is mostly u^, and standard z/^ ^ Uj. oscillations governed 
by Ai3 remain virtually undisturbed for atmospheric-like neutrino energies. 



III. CONSTRAINTS FROM THE NEUTRINO DATA 

Here we explore choices for the LIV parameters that may allow one to accommodate all 
neutrino data. We do not aim at performing a global fit to all data, but concentrate on 
qualitatively understanding and satisfying the relevant constraints. 

For LSND-hke energies and baselines {E ~ 50 MeV and L ~ 30 m), the v^. 
transition probability is easy to compute. At this point, we further assume that the energy 
dependence is such that, at least for LSND-hke energies, 

where a at is the smallest nonzero a„ term. The reason for this will become clear when we 
look at the solar neutrino data. We will comment on the behavior of the sum at large 
energies in due time. It is easy to compute all P^q in the limit |Ai3|L <^ 1: 

(6) 
(7) 

P^^ = 1 - 4 cos^ C sin^ Ol (1 - cos' C sin' Ol) sin' ( «iv (^-^ j - j • (8) 









1 2 


E^ 


2N-1 T - 




1 h 



In order to fit LSND and KARMEN-2 data, we must make sure that, at LSND, the 
oscillatory effects do not average out. Roughly speaking, we would like to choose a at so that 

Pe^^/Pf^^'^^^^ ~ (^LSND/^KARMEN)2^ ^ -g SUCCCSSful maSS-iuduCCd OSCiUatiOU 



6 



fits to LSND and KARMEN data [17\. Furthermore, we will assume that at experiments 
with larger energies and longer baselines, oscillation effects average out. Here, we will be 
particularl y co ncerned with^led searches for z/^ z/g and at NOMAD 18. loj]. 

CHORUS 20|, and NuTeV 121j. These constrain, at the 90% confidence level. 



sin^ 29 L cos^ C < 1-1 x 10-^ and sin^ 2C sin^ Ql < 3.3 x 10"^ 



(9) 



while the LSND data require sin^ 20 r, cos^ C ~ 10^'^. We therefore choose sin^2(^ = to 
avoid the NOMAD constraint 19], and are forced to choose sin^2^i = 1.1 x 10^'^ in order 
to avoid the NuTeV constraint |j2l| and to maintain hope that there is a passable fit to the 
LSND data. 

Solar neutrino data, combined with those from KamLAND, are very well explained by 
solar matter-affected v^. — > oscillations (where is some linear combination of and 
v-r)- This is clearly visible in the fact that, for ®B neutrinos. Pee is significantly less than 
one half, which, given our current understanding of neutrino masses and mixing, can only 
be obtained if non-negligible solar matter effects are at work. In order to avoid spoiling this 
picture, we require that, for solar-like neutrino energies, LIV effects are negligible. 

An order of magnitude estimate can be readily performed. The matter potential in the 
sun's core is of order A^an ~ 5 x 10^^ eV^/MeV, so that significantly smaller a^v values 
guarantee that, for large iV, Pee is undisturbed for energies less than Eq. This can also be 
deduce from Fig. ^ where the survival probability of solar electron neutrinos as a function 
of energy, for Eq = I'd MeV and = 5, is shown. The different curves correspond to 
05=0 ("normal" LMA behavior), ag = 5 x 10"^ eVVMeV, ag = 5 x 10"^ eVVMeV, and 
as = 5 X 10~^ eV^/MeV. The relevant standard neutrino parameters were set to Am^2 = 
8 X 10~^ eV^, sin^ 612 = 0.3, and sin^ 6*13 = 0. In order to compute Pee, we use the formalism 
outlined in j22], to which we refer for all the relevant details. 

Since the z/SM successfully describes solar neutrinos up to the "end" of the *B neutrino 
spectrum (around 15 MeV), its predictions should not be dramatically modified for these 
neutrino energies. Then, from the figure, we see that, for Eq = 15 MeV, solar data require 
05 ~ 10~® eV^/MeV. This estimate is qualitatively independent of A^. For larger values of 
A^, the same behavior is observed, except that the transition between LMA-like behavior 
and Pee = 1 is more abrupt. We conclude that if we choose Eq ~ 15 MeV and a^v ~ 
10~^ eV^/MeV, observable effects in solar data should be safely absent. Once such parameter 



2 4 6 8 10 12 14 16 18 20 

E (MeV) 



FIG. 1: Survival probability of solar neutrinos (Pee) as a function of the solar neutrino energy 
{E), for Eq = 15 MeV and different values of a^. The different curves correspond to = 5 and 
05 = (solid, black hue), 05 = 5 x IQ-^ eVVMeV (dashed, red hue), 05 = 5 x IQ-^ eVVMeV 
(dotted, green line), and 5 x 10^^ eV^/MeV (dashed-dotted, blue line). See text for details. 

choices are made, KamLAND data is completely oblivious to LIV effects, given that relevant 
reactor neutrino energies are less than 10 MeV. 

Next we consider atmospheric neutrinos. For the values of in which we are interested, 
LIV effects for E > 100 MeV are very large and ul is basically decoupled from the other 
two neutrinos. Thus, we are left with a standard effective two-flavor oscillation scheme 
between these two states, governed by the oscillation frequency A13. Since we choose ul to 
be mostly z/e, the other two neutrinos are mostly linear combinations of z/^ and z/,-, such that 
i^fi oscillations proceed as if there were no LIV effects. Effects due to the deviation 

of vl from a pure z/e state can be readily computed, and turn out to be very small. We 
have checked that, for the parameter values considered here, such deviations are within the 
current experimental uncertainties.^ 

All considerations above constrain 



so that A^ is the only left-over parameter. Since we wish to avoid averaged-out oscillations for 
typical LSND energies and baselines (if at all possible), the best one can do is to choose the 
oscillation phase to be close 7r/2 for typical LSND parameters, L ~ 30 m and ~ 50 MeV. 
This translates into 2N — 1 = 9, or A^ = 5.^ This is far from providing a good fit to the 

^ These effects are indeed very small, and it is not clear whether one will be able to probe them in future 

atmospheric neutrino experiments. 
^ For smaller values of N ^ there is "no hope" of reconciling solar and LSND data. A detailed discussion 




2N-1 



(10) 



8 



LSND data, for a couple of reasons. 

pLSND,max < IQ-S^ ^hich, even in the case of constant 
P^g^^°, already provides a mediocre explanation to the LSND anomaly. To add insult to 
injury, because of the acute energy dependence, Pj^^^^ falls very rapidly with energy, so 
that, on average, P^e ^ 10"^. 

Even if one overlooks the fact that we failed to accommodate all neutrino data, there 
is another issue we need to deal with. For values of the LIV parameters discussed above, 
and assuming that all other to vanish, the neutrino dispersion relation starts to differ 
significantly from that of a standard, ultra-relativistic particle for \p\ values larger than a 
few GeV. One way to avoid as uncomfortable a situation is to postulate that, for large values 
of IpI^, / either stops increasing or decreases. We will now argue that both possibilities not 
only avoid a dramatic modification of the neutrino dispersion relation for GeV energies, but 
allow one to contemplate accommodating all neutrino data more comfortably. 

We will briefly investigate two examples. One is to assume that 

where I is a positive integer, k is a constant coefficient and we chose to express the energy 
in "units" of 15 MeV, as before. In this case, the LIV potential is at most equal to k (in the 
limit E ^ 15 MeV), and we choose n so that, for LSND-like neutrino energies and baselines, 
kL ~ 1, which leads to k ~ 0.02 eV^/MeV. In order to accommodate the solar data, we need 
to choose / so that Ktanh^'^^(l) ^ 10"^ eV^/MeV. This translates into / ^ 16. It remains 
to discuss the faith of NuTeV (plus CHORUS and NOMAD) neutrinos. If we fix C = 0, P^r 
vanishes and, since here E ^ 15 MeV, 

K, \ f L 



P^e — sin^ 29 L sin^ 



5.1 



(12) 



0.02 eVVmeV ) VlOO m. 

i.e., the oscillation length is energy independent. In reality, we anticipate that oscillatory 
effects average out because different neutrinos propagate different distances between produc- 
tion and detection. Given that decay tunnels and detectors are several tens of meters long, 
and the oscillation lengths we are interested in are of order 30 m, it sounds like a reasonable 
assumption. Under this assumption the constraints listed in Eq. Q apply.^ 

was presented for TV = 1/2 (energy independent LIV effect) in [2^ . 
^ If oscillation effects did not average out, we would be allowed to choose k so that to Vf. oscillation 
effects at these experiments are minimized, hence loosening the constraints. 
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FIG. 2: Ue transition probability [P^e) as a function of energy {E) at (a) KARMEN (L = 

18 m) and (b) LSND (L = 30 m), using the tanh ansatz (Eg. compared with those obtained 
with marginal mass- induces oscillation solutions ((c) at KARMEN and (d) at LSND). See text for 
details. 

Fig. 121 depicts P^e as a function of energy for (a) L = 18 m (typical of KARMEN 
baselines) and (b) L = 30 m (typical of LSND baselines) for k = 0.02 eVVMeV, / = 16, 
and sin^ = 1.1 x 10^'^. Also depicted ((c)-|-(d)) are equivalent two-flavor mass-induced 
oscillations for Am^ = 1.5 eV^ and sin^26' = 1.5 x 10~'^, parameter choices that should 



fit the combined LSND+KARMEN data at the 99% confidence level {V^J This choice of 
parameters almost provides a good fit to all data, were it not for the fact that sin^ 29^ is 
constrained to be a little too small. 

A much more satisfactory fit could be obtained if the constraints listed in Eq. Q did 
not apply. In that case a good fit to the LSND data can already be obtained for sin^ 29i ^ 
2 X 10"^. It is plausible that a more detailed analysis of NuTeV, NOMAD, and CHORUS 
data, beyond the intentions of this short paper, would reveal that it is possible to loosen 
Eq. © by some factor of order one, but understanding whether this is indeed the case is far 
from trivial. Among the issues that need to be considered carefully is the fact that detectors 
and decay pipes are not negligibly small compared with the typical baselines, and one needs 
to average over their dimensions taking into account the exponential decay profile of the 
parent pions and kaons that originate the different neutrino beams. 

Finally, we mention the most radical change. Instead of considering LIV effects that 
behave like Eq. (fTT|) . we postulate that these are best represented by a function that is 



^ Recall, however, that the smallest allowed value of sin^ 29 at a given confidence level depends on the 
statistical analysis method employed [T^. 
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negligible at ^ 15 MeV and E ^ 1 GeV, while relevant for LSND energies "in between." 
In this case, almost by construction, all neutrino data can be comfortably accommodated. 
There is one concern that arises from atmospheric data. If |/| falls for energy values above 
100 MeV, at some energy one would encounter resonant to z/g transitions. While the 
resonance is expected to be very narrow and hence "harmless," it can be avoided altogether 
if one chooses ( and 6l so that ul is orthogonal to u^. 

Clearly, this possibility is even more tailor-made toward resolving the LSND anomaly. 
Moreover, as we will briefly comment later, we have no simple argument that justifies sup- 
pressing LIV effects at large energies other than very nontrivial fine-tuning. 



IV. SUMMARY AND REMARKS 



We have explored a new three-flavor fit to all neutrino data, including those from the 
LSND experiment. We postulate that one of the neutrinos couples to a source of Lorentz- 
invariance violation (LIV), so that its dispersion relation is modified in such a way that 
oscillation effects proportional to sin^(i?"L) are present for LSND neutrinos (see Eq. Q). 
For positive n values, we see that it is possible to accommodate solar and LSND data, while 
atmospheric data is easily accommodated as long as the LSND mixing angle and "oscillation 
length" are small. Curiously enough, the strongest constraints to this approach are provided 
by searches for z/^ Up oscillations at "small" L/E. Constraints obtained by NOMAD, 



CHORUS and NuTeV 



2l\ seem to allow only a mediocre fit to all neutrino data. 



unless the LIV parameters are chosen so that LIV oscillations are only relevant at LSND-like 
energies (30 MeV < E < 100 MeV). 

The model explored here is, by construction, very finely-tuned. Indeed, all manifestations 
of it are. The point where the LIV term becomes important is immediately below LSND 
energies and above typical solar neutrino energies. We were forced to choose very steep 
functions of the neutrino energy for energies between few MeV to tens of MeV. Clearly, this 
class of models is not motivated, but just serves as an example of how nontrivial it is to 
accommodate the LSND anomaly when the rest of the world neutrino data is included. In 
order to obtain a "natural" fit, we are all but forced into postulating that LIV effects are 
only present at LSND. 

The character of the LSND anomaly will, hopefully, be clarified with new experimental 
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data. The MiniBooNE experiment [23], currently taking data at Fermilab, was constructed 
in order to test whether the LSND anomaly is a manifestation of neutrino flavor transitions. 
As far as the different scenarios discussed here are concerned, MiniBooNE, which studies 
muon-type neutrinos with energies of several hundreds of MeV, is expected see either no 
"new physics" effect - if the LIV are indeed concentrated at LSND energies only - or a 
rather small signal consistent with "averaged out" oscillations. On the other hand, 

if MiniBooNE were to see a significant, energy dependent signal for z/g appearance, LIV 
models of the type discussed here would be significantly disfavored. Hence, in the case of no 
(significant) sign of new physics at MiniBooNE, we would be unable to rule out the kind of 
LIV effects discussed here. Conclusive information capable of resolving the LSND anomaly 
and addressing the models discussed here could be obtained, on the other hand, in new 
experiments with lower energy neutrinos, such as "long-baseline" studies of muon or pion 



decay at rest (see, for example, |2J]). 

We find it remarkable that so many distinct attempts to address the LSND anomaly 
by augmenting the z/SM fail because of distinct aspects of the world neutrino data (see, 
for example, jlJ], for a nice overview). Some L/E effects, for example, are disfavored by 
atmospheric and solar data (as in the case of "2+2" spectra) or by combined failed searches 
for and Ve disappearance at "short" baselines (as in the case of "3+1" spectra). L- 
independent solutions are disfavored by the KARMEN experiment (see, for example, \2l 



3V failed 



for a specific example), and we find that L x effects are ultimately disfavored 
searches at NOMAD, CHORUS, and NuTeV (this was also mentioned in 
A few model-building comments are in order. We need to introduce a term that violates 
Lorentz invariance and couples to neutrinos. Instead of providing a full model, we simply 
make the following general remarks. Consider a scalar field (j) whose time-derivative (in the 
preferred reference frame) acquires a vacuum expectation value: 



{d,<P) = 6o^M, (13) 

and we further impose a symmetry on 0. This symmetry may be softly broken in the 
potential, but is assumed to be (almost) exact when considering the couplings to neutrinos, 
so that only terms that scale like i?^" (and its powers) are allowed. Hence, symmetry 
arguments can be used to explain why only large values contribute. A similar argument 
applies for LIV effects proportional to [tanh(ii^)]'" as far as the exponent m is concerned, but 
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we cannot use symmetry arguments to "explain" why infinite series of liiglier dimensional 
LIV operators would combine into either a tanh(i?) form, or why effects should be enhanced 
at neutrino energies close to 50 MeV (but severely suppressed at energies above or below 
50 MeV). 

In summary, we succeed in identifying a three-neutrino, LIV solution capable of accom- 
modating all neutrino data. The price we had to pay in order to do that was quite high: not 
only were we forced to impose that the LIV effects are very strongly energy-dependent for 
energies around 20 MeV, we were also forced to postulate that LIV effects either "flatten 
out" for energies larger than 100 MeV (in which case the fit to LSND data is mediocre) or 
peak at LSND energies and quickly "disappear" ioi E ^ 1 GeV. 

We do not advocate our model as a solution to the LSND anomaly, but rather as another 
evidence for the "real" LSND puzzle: how can the effect observed at LSND be due to physics 
beyond the z/SM if all other experiments ("neutrino" or otherwise) have failed to observe 
related "new physics" effects? If the integrated literature on the subject is to be taken as 
a good indicator, one would, naively, conclude that the answer to the question above is 
'it cannot.' It seems unfair, however, to conclude that a new physics interpretation of the 
LSND anomaly does not exist simply because there seems to be no elegant solution to the 
LSND anomaly. It may well be we just have not found one yet. 
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